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In classical general relativity the Cauchy horizon within a two-horizon black hole is unstable via a 
phenomenon known as mass inflation, in which the mass parameter (and the spacetime curvature) 
of the black hole diverges at the Cauchy horizon. Here we study this effect for loop black holes 
- quantum gravitationally corrected black holes from loop quantum gravity - whose construction 
alleviates the r = singularity present in their classical counterparts. We use a simplified model of 
mass inflation, which makes use of the generalized DTR relation, to conclude that the Cauchy horizon 
of loop black holes indeed results in a curvature singularity similar to that found in classical black 
holes. The DTR relation is of particular utility in the loop black hole because it does not directly rely 
upon Einstein's field equations. We elucidate some of the interesting and counterintuitive properties 
of the loop black hole, and corroborate our results using an alternate model of mass inflation due 
to Ori. 
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I. INTRODUCTION 

A significant hinderance to the development of the cor- 
rect theory of quantum gravity lies in the fact that, by 
its very nature, any predicted phenomenology is likely 
to be beyond the reach of practical experimental capa- 
bility. Without any experimental guidance (or at least 
very little) it seems that we are merely left to postulate 
how quantum gravity may alter classical theory in high 
energy regimes, with constraints on the theory coming in 
the form of self-consistency checks. Despite this limita- 
tion there are still predictions that can be used to guide 



us in our quest to find the fundamental theory. Of partic- 
ular interest in this respect is how quantum gravity may 
change the presence of the singularities that are predicted 
by general relativity. It has long been expected that the 
infinities that occur, for example inside black holes, rep- 
resent a breakdown of Einstein's theory and that a proper 
quantum gravitational treatment will be able to placate 
them. This then provides us with a "laboratory" that 
can be used to test the merit of a proposed theory of 
quantum gravity; for example a theory predicting that 
the interiors of black holes are regular can be regarded 
as an improvement over those that do not. 

With this as motivation, much work has been per- 
formed to say how black holes might be modified in the 
context of quantum gravity. Several proposals from sev- 
eral different approaches to quantum gravity have been 
given that are free of the classical r = singularity [THS] • 
Despite being regular at r = however, these quantum 
black holes generically have more than one horizon. One 
therefore ought to be wary of their overall regularity since 
it is well known that classical black holes with more than 
one horizon also exhibit divergent behaviour at their in- 
ner horizon (the Cauchy horizon) . A demonstration that 
quantum black holes also exhibit divergent behaviour at 
their Cauchy horizon will indicate that they are not per- 
forming the job that is expected of them; i.e. they do not 
represent regular black holes. At the least this will imply 
some insufficiency in the assumptions and (semiclassical) 
approximations used in constructing the solutions, if not 
in the underlying theory of quantum gravity from whence 
they came. 

In classical gravity, the instability of the Cauchy hori- 
zon can be seen to arise from the infinite blueshift of 
external radiation that occurs on the inner horizon. This 
phenomenon is most easily demonstrated in the Reissner- 
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Nordstrom black hole, but can also be shown to arise 
in the Kerr and Kerr-Newmann black holes. This type 
of instability had been thought to occur for many years 
and was demonstrated to first order by observing the 
behaviour of perturbations on the metric (HJ [7J. Such 
an analysis, however, says nothing about the backreac- 
tion on the spacetime that is induced. It was not until 
1990 that Poisson and Israel showed rigorously that in 
a realistic black hole the Cauchy horizon of a Reissner- 
Nordstrom spacetime is unstable via a phenomenon that 
they dubbed mass inflation. The Cauchy horizon generi- 
cally develops singular curvature caused by the local mass 
parameter that grows indefinitely large [S]- 

In this paper we ask whether or not mass inflation 
occurs in a class of quantum gravitational black holes 
motivated by loop quantum gravity. We will see that 
a mass-inflation-like phenomenon indeed seems to oc- 
cur. Although we work within the context of a specific 
model, our results depend on quite general features of the 
black hole interior, suggesting that any black hole with 
a Cauchy horizon - classical or quantum - will be hard 
pressed to avoid the same result. 

Loop quantum gravity (LQG) (TUlU2| . has given rise 
to models that afford a description of the very early uni- 
verse. This simplified framework, which uses a minisu- 
perspace approximation, has been shown to resolve the 
initial singularity problem |13) . A black hole metric in 
this model, known as the loop black hole (LBH) [T], was 
constructed from a modification to the holonomic version 
of the Hamiltonian constraint. A parametric function 
<j(6), that labels the elements in the class of Hamiltonian 
constraints compatible with spherical symmetry and ho- 
mogeneity, can be uniquely fixed from asymptotic bound- 
ary conditions, yielding the proper classical Hamiltonian 
in the limit the polymeric parameter <5 — > 0. 

As such, this black hole solution might be expected to 
model features induced by LQG effects. It has a property 
of self-duality that removes the r = singularity and re- 
places it with another asymptotically flat region. Both 
the thermodynamic properties [TJ 114) and the dynami- 
cal aspects of collapse and evaporation [T5] of these self- 
dual black holes have been previously studied. However, 
since the polymerization of the Hamiltonian constraint 
in the homogeneous region is inside the event horizon, 
the physical meaning of the solution when the metric is 
analytically continued to all spacetime remains an open 
problem. Nevertheless this metric is useful insofar as it 
can be expected to provide a first approximation to black 
hole solutions that emerge from LQG. Since its interior 
has a Cauchy horizon, we can investigate its stability in 
this quantum gravitational framework. Black hole space- 
times have also been investigated in a midi-superspace 
reduction of LQG |T7j . While the black hole solution ob- 
tained by this method is closer to the full LQG theory, 
it can only be presented in numerical form, whereas the 
black hole we consider has a closed analytic form and so 
can be more easily investigated. We emphasize that our 



analysis has limitations insofar as it does not correspond 
to a full LQG solution. 

Our paper is organized as follows. In Sect, [n] we re- 
call the loop black hole (LBH) derivation in short and 
expound some of its properties. In Sect. |HT] we con- 
sider a short perturbative calculation that indicates the 
instability of the Cauchy horizon and motivates further 
analysis. In Sect. |IV| we will review the derivation of 
the so called generalized DTR relation. This result will 
play a vital role in Sect. [V| where we will see how mass 
inflation arises in the loop black hole before concluding 
in Sect. |VI| We consider in an Appendix a model of 
mass inflation developed by Ori [TB] applied to the LBH. 
Notwithstanding certain difficulties entailed in this ap- 
proach, we find that it corroborates the main results in 
our paper. 



II. THE LOOP BLACK HOLE 

The regular black hole metric that we will be using 
is derived from a simplified model of LQG pQ. LQG is 
based on a canonical quantization of the Einstein equa- 
tions written in terms of the Ashtekar variables [TS] , that 
is in terms of an SU(2) 3-dimcnsional connection A and a 
triad E. The basis states of LQG then are closed graphs, 
the edges of which are labeled by irreducible SU(2) rep- 
resentations and the vertices by SU(2) intertwiners (for 
a review see e.g. [TUHT2"] ). The edges of the graph repre- 
sent quanta of area with area 7^p\/j(j + 1), where j is 
a half- integer representation label on the edge, lp is the 
Planck length, and 7 is a parameter of order 1 called the 
Immirzi parameter. The vertices of the graph represent 
quanta of 3-volume. One important consequence that 
we will use in the following is that the area is quantized 
and the smallest possible quanta correspond to an area 
of V3/2jl 2 p . 

To obtain the simplified black hole model the following 
assumptions were made. First, the number of variables 
was reduced by assuming spherical symmetry. Then, 
instead of all possible closed graphs, a regular lattice 
with edge lengths 5b and S c was used. The solution was 
then obtained dynamically inside the homogeneous re- 
gion (that is inside the horizon where space is homoge- 
neous but not static, a spacetime known as Kantowski- 
Sachs spacetime). An analytic continuation to the out- 
side of the horizon shows that one can reduce the two 
free parameters by identifying the minimum area present 
in the solution with the minimum area of LQG. The one 
remaining unknown constant Sf, is a dimensionless pa- 
rameter of the model that determines the strength of 
deviations from the classical theory, and would have to 
be constrained by experiment. Redefining Sf, — <5, the 
free parameter that appears in the metric is e = ^7, 
the product of the Immirzi parameter 7 and the poly- 
meric quantity S. With the plausible expectation that 
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quantum gravitational corrections become relevant only 
when the curvature is in the Planckian regime, corre- 
sponding to ^7 < 1, outside the horizon the solution is 
the Schwarzschild solution up to negligible Planck-scale 
corrections in lp and Sj. 

The procedure to obtain the metric is the following. 

i We define the Hamiltonian constraint by replacing 
the homogeneous connection with the holonomies 
along the fixed graph decided above. The diff- 
constraint is identically zero because of homo- 
geneity and the Gauss constraint is zero for the 
Kantowski-Sachs spacetime. 

ii We solve the Hamilton equation of motion for 
the holonomic Hamiltonian system imposing the 
Hamiltonian constraint to be zero. 

iii The third step consists in extending the solution to 
all spacetime. Note that while there is no math- 
ematical obstruction to this extension it is an as- 
sumption, since we cannot know the correct poly- 
merization of the Hamiltonian constraint in the full 
spacetime as we found the solution in the homoge- 
neous region. 

A more technical point is the following. In the homo- 
geneous region a cut-off in the spatial direction is intro- 
duced. This extra structure disappears in the metric, 
and all physical quantities are independent of this length 
scale pp. 

The above procedure then yields a black hole solution 
that can reasonably be expected to model features in- 
duced by LQG effects. This quantum gravitationally cor- 
rected Schwarzschild metric - known as the Loop Black 
Hole (LBH) - can be expressed in the form 



ds 2 = -G{r)dt 2 
with dVt 2 = d6 2 + sin 2 
G(r) 



dr 2 
2 and 



H(r)dn 2 



(1) 



(r - r + )(r - r_)(r + r. t )' 2 



F(r) 



H(r) 



(r — r + )(r — r_)r 4 
(r + r,) 2 (r 4 +a2) 

„2 



(2) 



Here, r + = 2m and r_ = 2mP 2 are the two horizons, 
and = ^/ryrZ — 2mP. P is the polymeric parameter 

P = (\/l + e 2 -l)/(\/l + e 2 +l), with e = <fy < 1. Hence 
P -C 1, implying r_ and are very close to r = 0. The 
area a is equal to j4 m ; n /87r, A mm being the minimum 
area gap of LQG. 

Note that in the above metric, r is only asymptotically 
the usual radial coordinate since ggg is not just r 2 . This 



choice of coordinates has the advantage of easily reveal- 
ing the properties of this metric. Most importantly, in 
the limit r — > oo, the deviations from the Schwarzschild- 
solution are of order Me 2 /r, where M is the usual ADM- 
mass: 



G(r) 

F(r) 
H(r) 



2M , , 

1 (1 - e 2 

r 



1 - 

,,.2 



2M 
r 



(3) 



The ADM mass is the mass inferred by an observer at flat 
asymptotic infinity; it is determined solely by the metric 
at asymptotic infinity. The parameter m in the solution 
is related to the mass M by M = m(l + P) 2 . 

Performing the coordinate transformation R = a /r 
with the rescaling t = tr 2 /a a , and simultaneous replace- 
ments R± = a /r T , R* = a /r* yields the result that 
the metric in the new coordinates has the same form as 
that in the old coordinates, thus exhibiting a very com- 
pelling type of self-duality with dual radius r = ■ s fK . 
Looking at the angular part of the metric, one sees that 
this dual radius corresponds to a minimal possible sur- 
face element. It is then also clear that in the limit r — > 0, 
corresponding to R — > oo, the solution does not have a 
singularity, but instead has another asymptotically flat 
Schwarzschild region. 

Computing the surface gravity 

k 2 = W„x72 = -W%LW2, (4) 

where — (L 0, 0, 0) is a timelike Killing vector in r > 
r+ and r < r_ (but spacelike in r_ < r < r + ) yields 



4m 3 P 4 (l - P 2 ) 



4m 3 (l - P 2 ) 
16m 4 + a 2 



(5) 



for the surface gravity on the inner and outer horizons of 
the metric 0. 

We now go on to a simple calculation that both in- 
troduces some of the background structure necessary for 
investigating mass inflation and gives a first illustration 
of the troublesome behaviour exhibited by the Cauchy 
horizon. 



III. A PRELIMINARY CALCULATION OF 
STABILITY 

In this section we are going to perform a 'warm up' 
calculation for mass inflation. We start with the LBH 
given by Eq. ([!]), though the same calculation is eas- 
ily performed for any spherically symmetric black hole. 
Consider a stream of null dust entering the black hole. 
We neglect the backreaction that this radiation has on 
the spacetime and we compute the energy density that 
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a timelike observer would measure from it as they cross 
the Cauchy horizon. 

We will find it useful to define the advanced time 
v = t + r*, where r* is the tortoise coordinate (not to 
be confused with the variable = 2mP present in the 
metric) defined by 



dr* 
dr 



1 



(6) 



The limit v = oo corresponds to future null infinity when 
outside of the black hole and to the Cauchy horizon when 
inside. Thus, when we speak of "approaching the Cauchy 
horizon" we are referring to the v — ¥ oo limit. It has 
been shown [TI5] that for a black hole formed by gravi- 
tational collapse there will always be ingoing null radia- 
tion with energy density falling off as an inverse power, 
~ v^ 1 where 7 > 12, for large v. In any gravitational 
collapse there will always be gravitational radiation that 
escapes into space, but some of this radiation will later 
be scattered back by the gravitational potential present 
in the spacetime surrounding the black hole. The inverse 
power law given here corresponds to the luminosity of 
this backscattered radiation, and its importance lies in 
the fact that it is expected to always be present in a real- 
istic black hole. We make use of this result for the LBH, 
since the gravitational emission and backscatter in ques- 
tion occurs outside of the event horizon, where the LBH 
spacetime only negligibly deviates from that of a classi- 
cal Schwarzschild black hole, for which Price's result is 
known to apply. This situation is depicted in Fig. [T] 

First we cast the loop metric in (v, r, 9, </>) coordinates, 
obtaining 



ds 2 = -G{r)dv 2 




drdv + H(r)dn 2 , (7) 



where the metric functions are as defined in Sect. ITT1 
We then consider a stress-energy tensor corresponding 
to ingoing null dust. It takes the form 



T a fj = fJ,(r,v)(d a v)(dpv). 



(8) 



With these last two equations we determine the form 
of /x(r, v) by computing the conservation condition = 
T a fn. The result is that we must have /i(Y, v) = A(r)L(v), 
where L{v) is an arbitrary function of v and A(r) must 
satisfy the equation 



= r(r 4 + a 2 ) dA{r) > ^ 



dr 



2(r*-a 2 )A(r), 



the solution of which is 



A(r) oc 




FIG. 1: A Penrose diagram displaying the spacetime 

outside of the black hole as well as between the 
horizons. We illustrate the ingoing null radiation and 
label the direction of increasing advanced time v. 



and so we find 



fj.(r,v) 



4ir(r 4 



-L(v) 



(9) 



for the general form of /x(r, v). 

Far from the black hole the spacetime is well ap- 
proximated by a Schwarzschild solution, and indeed for 
r » y/ao this form of /i(r, v) reverts to the same form as 
is known for the Schwarzschild black hole. We thus iden- 
tify L(v) with the usual luminosity, and we expect the 
analysis by Price [T!5] to hold to a good approximation. 
That is, we have L(v) ~ w~ 7 with 7 > 12. 

In order to determine the energy density a timelike 
observer would measure from this dust we also need to 
specify the four-velocity of the observer. For the case at 
hand, however, all we need is the v component, which we 
find to be 



JJV 



(E-y/EP-G), 



(10) 



where E is a constant of the motion, the sign of which 
gives the direction of travel between the horizons. E > 
corresponds to an ingoing (left-moving) observer and 
E < to an outgoing one (right-moving). We want to 
focus on right-moving observers, since it is they who will 
reach the right branch of the Cauchy horizon (the branch 
located at v = 00); we thus only consider geodesies with 
negative E. The energy density measured by such an 
observer is given by 



p = T a pU a UP =fi( ri v)(U v ) 2 



(11) 
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As the Cauchy horizon is approached we have G — > 0, 
and so the negativity of E implies that 



{U v f ~ 4S 2 /G 2 , 



(12) 



as this approach takes place. However, it is clear that we 
also need to know in what manner the metric function 
approaches zero. For this one can compute the tortoise 
coordinate r* directly from its definition and show that 
near the Cauchy horizon, r ~ r_, it is dominated by a 
term given by 

r* ~ -(2K_)- 1 log|r-r_|, 

where k_ is the surface gravity of the Cauchy horizon, 
given by Eq. (fjj). If one then computes the t and r 
components of the four- velocity it can be shown that for a 
right-moving observer = ^ ~ 1 near the Cauchy 
horizon. From this we obtain 



— v = —t — r* ~ — 2r* + const 

~ kZ 1 log |r — r_ | + const => (r — r_) _1 oc e 14 "". 

which, since the metric functions F and G are propor- 
tional to r — r_ , results in the conclusion that near the 
Cauchy horizon they decay as 



G and F oc e 



oo. 



(13) 



Thus, we finally find that the energy density measured 
by the observer goes as 



9 1 



- V L( 



(14) 



as the Cauchy horizon is approached. Since L(v) fol- 
lows an inverse power law, we conclude that the observer 
measures diverging energy density as they cross of the 
Cauchy horizon. 

A simpler scenario, in which a field pulse is fired in 
to the black hole and then scattered in the interior ge- 
ometry, yields a different result [8|. Its energy density 
diverges at the Cauchy horizon in a manner dependent 
on the difference of the surface gravities k_ and k + . Pro- 
vided k + > K- , no divergence occurs, and the loop black 
hole appears to be more stable under this type of radia- 
tion than classical black holes. However there does not 
appear to be any cure (within the perturbative frame- 
work) for the more generic result obtained above using 
inverse power law decay. 

With this result as motivation we move on to the mass 
inflation analysis, computing the behaviour of the space- 
time curvature resulting from the infinitely blue-shifted 
radiation. We employ a simplified but well-motivated 
model that makes the problem tractable. In preparation 
for this we will briefly review the derivation of the gen- 
eralized DTR relation, a key ingredient in our analysis. 



IV. THE GENERALIZED DTR RELATION 

One of the key ingredients needed for mass inflation to 
occur is the existence of outgoing null radiation within 
the black hole, in addition to the ingoing radiation (that 
which is undergoing infinite blueshift). It is the cross flow 
between these two streams that induces the phenomenon. 
It is expected that in a realistic black hole (one formed 
by gravitational collapse) such outgoing radiation will 
always be present, having originated from the surface of 
the collapsing star. Performing this calculation for the 
LBH using continuous streams of ingoing and outgoing 
radiation as done originally (9 appears to be technically 
intractable. Fortunately we can make use of a model sys- 
tem that contains the same essential physics but is much 
easier to deal with technically. Specifically, we replace 
the continuous ingoing and outgoing streams by null thin 
shells that collide near the Cauchy horizon. Such an anal- 
ysis was performed, for example, to show that mass infla- 
tion occurs in rotating black holes [2U] (though the rotat- 
ing case was first considered using the Ori model [H]). A 
further motivation for using this simplified model is that 
it allows us to make conclusions about the spacetime cur- 
vature without having to directly use the Einstein field 
equations. In fact the only reference to field equations 
that we will make is to assume that the Einstein equa- 
tions hold to good accuracy far from the black hole; this is 
of course a well justified assumption. Not having to rely 
directly on field equations affords a significant advantage, 
because we wish to test stability given a loop quantum 
gravity framework, a framework in which we have limited 
knowledge of the effective field equations that should be 
used. Of course this analysis is still semiclassical, but 
we should try to avoid using the classical field equations 
whenever possible, and the null shell model we will use 
here does just that. 

Colliding null shells act to split the spacetime into four 
separate regions, as illustrated in Fig. [2]. The DTR rela- 
tion is a result that allows us to relate the metrics in each 
of these regions to each other (evaluated on the collision 
two-sphere S). The original DTR relation achieves this 
but is constrained to Schwarzschild spacetimes only [22] . 
The general DTR relation [201 [23] , whose derivation we 
briefly recapitulate, evades this limitation. 

While we assume a spherically symmetric spacetime 
and shells that travel radially, the results presented are 
independent of this latter assumption. For a given radial 
null shell, E, we define the vector tangent to its genera- 
tors 



t a = 



dx a 



(15) 



where A parameterizes the shell; this parameter must be 
the same on both sides of E but is otherwise arbitrary 
and need not be afflne on either side (in fact if the shell 
in question exerts transverse pressure then A cannot be 
made affine on both sides [13]). We also define the two 
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FIG. 2: Two colliding null shells separate spacetime 
into four distinct regions: A, B, C and D. In general 
the shells before the collision, 3 and 4, will not have the 
same properties as the shells after collision, 1 and 2. 
The two-sphere on which they collide is labeled S. 



tangent vectors e" that are transverse to the generators 
(latin indices run over two coordinates). Since S is null 
we have 



(16) 



With our assumption of spherical symmetry the trans- 
verse vectors must lie in the two-sphere, and so can be 
chosen to be e^d a = dg and e^d a = without loss 
of generality. We can then define the effectively two- 
dimensional intrinsic metric on E by 



(17) 



The intrinsic metric for null shells is effectively two- 
dimensional because including the third tangent vector 
£ a gives trivially zero entries: <j\\ — a\ a — as can be 
seen from Eq. (161; we therefore exclude these entries 



from the definition of the intrinsic metric. 
We now define an object associated with E, 



K„ 



1 



-Sff {(Tab 



1 



%i (fa/3) e"e 



(18) 



where is the Lie derivative along the vector l a . We 
have dropped additional terms on the right-most side of 
the above equation since they evaluate to zero anyway, 
as a result of e° being transverse to A. We call the trace 
of this K, and it is equal to the expansion rate of E; we 
find 



K 



Kn 



K, 



(19) 



where k is the "acceleration" of the shell defined by 
nl a . 

In order to consider what happens on the intersec- 
tion two-sphere S we need to first make the principal 
assumption that the spacetime on S is no worse behaved 



than it is on the thin-shells individually. To understand 
what this means one must realize that by introducing 
an infinitely thin shell with nonzero stress-energy into 
our spacetime we are necessarily introducing a curvature 
singularity on the shell. This singularity simply arises 
as an artifact of the model and has a well understood 
physical interpretation. Despite this singularity, given 
a local coordinate system that is continuous across the 
shell the metric g a p remains continuous and piecewise 
continuously differentiable across the shell (the reader is 
referred to [231 [M] for expositions on the thin-shell for- 
malism). Thus, the assumption that the spacetime on S 
is no worse than it is on the shells individually is equiva- 
lent to requiring that for every point on S there is a local 
coordinate system in which the metric is continuous and 
piecewise continuously differentiable. 

This assumption leads to a well defined notion of par- 
allelism in vectors transverse to S. In our model we have 
four null vectors that are orthogonal to S: if = dx a /dXi, 
i = 1,2,3,4, where Aj parameterizes shell £j. And since 
there are only two linearly independent null directions 
orthogonal to S the above assumption allows us to con- 
clude that 



h = 



h = b£ 4 , 



for some numbers a and b. From this we immediately 
obtain 

(4^2)(4-4) = (^i-4)(^-4) on S (20) 

This is a necessary geometric condition on the spacetime 
and is independent of any field equation. 

Let us now define four functions on S, Za,...,Zd as, 
for example 



Z A 



K X K 2 



(21) 



where Ki is the expansion of shell S», as defined in Eq. 
(19). The other functions are similarly defined, where 



the subscript on Z corresponds to the spacetime region 
bounded by the two shells referred to on the right-hand- 
side. Of principle importance is the fact that these func- 
tions are independent of the parameters Aj used to de- 
scribe the shells. This is easily seen because under repa- 
rameterization Aj — > Aj the expansion Ki rescales the 
same as ti. This allows us to put Eq. (20) into a 
parameter-independent form 



\Z A Z 



A^B\ 



\Zr,Z 



C^D\ 



on S , 



(22) 



where the absolute value bars are necessary because in 
using the Z functions we have lost sign information. For 
example, Za is invariant under the reparameterization 
Ai — > — Ai, unlike t\. 

This is what is known as the generalized DTR rela- 
tion. It allows us to relate the metric in the four space- 
time quadrants (on S) in a way that is independent of 
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the parameters Xi used to generate the shells and indeed 
is completely independent of the stress-energy properties 
of the shells. That is, nowhere in the derivation of this 
result have we used Einstein's field equations; rather the 
relation is a necessary geometric condition of the space- 
time. 



V. MASS INFLATION 

Classically, mass inflation is a process in which the 
mass parameter of the black hole diverges unboundedly 
at the Cauchy horizon. The outgoing radiation acts to 
separate the inner apparent horizon of the black hole 
from its Cauchy horizon, and this separation is a nec- 
essary condition for mass inflation to occur. The effect 
is induced independent of whatever form the outgoing 
radiation may take; the only requirement is that it has 
a nonzero energy (necessarily negative energy since it is 
outgoing). That is, we expect the radiation to change 
the mass parameter of the black hole, thus creating the 
separation between the apparent and Cauchy horizons 
stated above. This process was originally explicated for 
Reissner-Nordstrom black holes [9] and then simplified 
in a model constructed by Ori [TB] (we will discuss this 
model of mass inflation in the Appendix). From a more 
physical viewpoint, the inability of the counter streams of 
radiation to locally travel faster than the speed of light 
relative to each other can be seen to produce arbitrar- 
ily large amounts of gravitational energy as the Cauchy 
horizon is approached, and this energy manifests itself in 
the material form of arbitrarily large local mass. Hamil- 
ton and Avelino provide an excellent discussion of the 
precise physics at work here [35] and they note that the 
simplified models of mass inflation typically used (includ- 
ing the one we employ) do well at accurately capturing 
the results obtained in their more comprehensive study. 

The calculation that we will perform models the ingo- 
ing and outgoing radiation as null thin shells that collide 
near the Cauchy horizon, as displayed in Fig. [3] This 
model captures the same essential physics present in the 
continuous scenario while making the problem much sim- 
pler to solve. Furthermore, it allows us to make conclu- 
sions without the need to rely on specific field equations, 
the form of which are not known given a loop quantum 
gravitational framework. The colliding shells act to split 
the spacetime into four distinct regions, each equipped 
with their own metric. A principal assumption we make 
is that each of the spacetime quadrants is described by 
a loop metric of the form given by Eq. |T]), each with a 
different mass parameter m. This assumption seems rea- 
sonable considering that the loop metric was derived as 
a modified spherically symmetric vacuum solution, and 
each of the four regions are themselves representative of 
a spherically symmetric LBH vacuum spacetime. It is 
also supported by the requirement that outside the black 




FIG. 3: The colliding null shells split the loop 
spacetime into four regions, the metrics in each of which 
can be related via the generalized DTR relation. We 
consider the scenario in which the ingoing shell is 
arbitrarily close to the Cauchy horizon. 



hole classical general relativity holds to a good approx- 
imation; in a Schwarzschild black hole endowed with an 
ingoing thin-shell we know that the mass parameter is 
constant within the past and future regions of the shell. 
Thus m in the loop black hole can also be expected to 
have this behaviour for large r. Physically we expect 
stationarity to hold, meaning that m should be indepen- 
dent of the advanced time v within this region. Hence in 
each quadrant there will be only small deviations from 
constant m. 

With this assumption we can use the generalized DTR 
relation, Eqs. (21 1 and (22), to say how the mass pa- 



rameters in each region, ro^, mg etc, are related to each 
other on the collision two-sphere S. Since the DTR rela- 
tion itself only gives us information on S, this assumption 
allows us to say that the parameter tub, for example, 
takes the same value on S as it does in the rest of the B 
region up to negligibly small corrections. 

We wish to find how the relation ( 22 ) evaluates given 



that each region is described by the loop metric. Working 
in (t, r, <j)) coordinates the tangent vectors e" transverse 
to the generators of the radially traveling shells are given 



by e£ = (0,0,1,0) and 



(0,0,0,1). The intrinsic 



metric on the shells, Eq. (17), is therefore given by 



a ab = H(r) 



1 
sin 2 6» 



H{r) = r 2 



(23) 



where ao is interpreted as the minimum area of the the- 
ory. Recall that i a is the vector tangent to the generators 
of the shell in question. The entries of the intrinsic met- 
ric are scalars with respect to the full four dimensional 
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spacetime, and so when we compute the Lie derivative of 
them along the vector l a we simply use a partial deriva- 
tive: 



Note that we have assumed continuous r at S: this 



K ab = (a ab ) = h a {<j ab ) a 



(24) 



From this it is easily seen that the trace, Eq. (19), takes 
the form 

K = -^£ a H, a , 

or rather, when referring to one of the specific shells Si 
we have 



K — —P a H 



(25) 



We now obtain the form of the quantities Za, Zg etc, as 



defined by Eq. (21 ). For example we have 

i e?4H ta H tP 



Za 



H 2 li ■ l-x 



(26) 



where £" and f% are the vectors tangent to the generators 
of shells Si and S2 that bound spacetime region A. 

To relate Za to the metric g"^ in region A we make 
use of the completeness relation 



a _ Ml *2 , ab a 8 
9(A)- ZaH- 



(27) 



This can be verified by computing all inner products be- 
tween I" , £2 and e" . If we then contract the metric with 
H^ a Hp the result is seen to be 



afi _ 2H KiK 2 

1 (A) H a H p - — . 



(28) 



where we have used e"H a = 0, since that H a = for 
all a 7^ r. Comparing this result with Eq.(26) we imme- 
diately obtain 

Z * = ^p9^H, a H, = ^g^iHrf. (29) 

For the loop black hole we have H = r 2 + a^/r 2 , and so 



_ 2 {r'-a 2 ) 2 



r 2 (r 4 + ag) 



(30) 



where Fa is the metric function g rr in region A. Under 
the assumption that the region is described by the loop 
metric this function is given by 



(r — 2m A )(r — 2m.AP 2 )r 4 
A = {r + 2m A P) 2 (r i + a 2 ) ' 



(31) 



where ttia is the mass parameter in region A. The forms 
of Z in the other regions follow analogously to this. 



results from our assumption in Sect. IV that the space- 
time at S is no worse behaved than it is across the shells. 
We know that r must be continuous across the shells in 
order for them to have a well defined intrinsic geometry 
(specifically, a ab must be continuous across the shells) 
and so the assumption that r is continuous on S follows 
directly from the assumption already made in deriving 
the generalized DTR relation. 

Recall that the generalized DTR relation tells us that 
the four regions of spacetime, as seen in Fig. [3|, are 
related at S by \Z A Z B \ = \Z C Z D \. From Eq. @, the 
form of the metric function F and the assumption of 
continuous r we finally obtain the defining equation of 
this analysis: 



\X A X 



B 



\X c: X 



where for convenience we have defined 



2m A )(r- 2m A P 2 ) 



(r + 2m A P) 2 



(32) 



(33) 



and similarly for Xg etc. This equation gives us a nec- 
essary condition on the four mass parameters 

We now wish to ask the question of what happens when 
the ingoing shell is arbitrarily close to the Cauchy hori- 
zon (equivalently, the shell enters the black hole arbi- 
trarily far in the future as seen by an outside observer). 
In order to do this we must specify how the energy of 
the shell should depend on its distance to the Cauchy 
horizon or, more precisely, how it should depend on the 
advanced time v. In Sect. Ill we discussed how, in the 



case of continuous radiation, the energy density of the in- 
going radiation should fall off as an inverse power, ~ 
where 7 > 12, in the large v limit (recall that the Cauchy 
horizon is located at v = 00). This occurs as a result of 
backscattered gravitational radiation, and is considered 
a very generic form of matter present in realistic black 
hole spacetimes. In classical spherically symmetric black 
holes this ingoing radiation acts to increase the mass of 
the black hole as a function of v such that the change in 
mass is proportional to the power law decrease expressed 
by the energy density; i.e. m'(v) oc v^ 1 . This classical 
result relies on Einstein's field equations, and so cannot 
necessarily be expected to hold in the loop black hole. 
However, it is known that far from the loop black hole 
the spacetime is negligibly deviant from that of a clas- 
sical Schwarzschild black hole and that Einstein's equa- 
tions are expected to hold to a good approximation in 
this region. Therefore, under the assumption that the 
mass parameter m of the loop metric is fixed along a line 
of constant v (in the null shell model this corresponds to 
assuming fixed m within a given quadrant) it means that 
the classical result of m'(v) oc v~ 7 will hold to a good 
approximation in the loop black hole as well. When it 
comes to the null shell model, however, we must still 
make a choice about how to model the energy of the in- 
going shell as a function of v. This leads us to our third 
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primary assumption: that the energy of the ingoing shell, 
and therefore the jump in mass parameters across the 
shell, follows this same power law 



771(7 ~ m B v 



(34) 



This is the most natural condition that can be placed 
on the shell energy given the nature of the model being 
considered, and it is what will most accurately capture 
the important physics of the full system. The reader 
should note that this is the one and only place that any 
reference to the Einstein held equations is made, and it is 
simply the statement that the field equations should hold 
when far from the black hole, which indeed they should. 

There is a source of possible confusion that should be 
addressed. Previously we stated that the mass parame- 
ters should be fixed within their given regions, and yet 
Eq. (34 1 seems to imply that they are changing. What 
must be understood is that the action of "moving the 
ingoing shell closer to the Cauchy horizon" really just 
means that we are thinking of a series of distinct space- 
times. A given spacetime will have the ingoing shell lo- 
cated at some fixed value of v, and in our model the 
difference in mass parameters mc — m B depends on this 
value via Eq. ( 34 ) . With this framework we are going to 



consider what happens when the spacetime in question 
has its shell located at an arbitrarily large value of v, but 
the mass parameters of this spacetime are still constant 
within their respective quadrants. For the rest of this 
section when something is said to be a function of v it is 
implied that v is the value of advanced time at which the 
shell is located, and by increasing v we moving through 
a sequence of spacetimes. The only way in which the 
masses can be said to change as a function of v is within 
this context. 

With this understood, it must be specified that mg, 
the mass of the black hole as detected by an outside ob- 
server prior to the shell entering the hole, is to be con- 
sidered "fixed" as we cycle through different spacetimes 
(as the shell's value of v is increased) and the mass after 
the shell has entered, mc, will "change" over the series 
of spacetimes. Specifically, as a function of v the future 
mass is seen to decrease as mc — mB +v~ 7 . This choice 
is both the most natural (as opposed, for example, to 
considering mc as fixed and ms as increasing) and it 
allows us to use a well defined choice of advanced time 
v. That is, the advanced time we use will be that which 
is defined with respect to the B region of the spacetime. 
Equivalently, the choice of fixed m B tells us that that 
the Cauchy horizon will be located exactly at r = r_ 
and v — oo, where 



2m B P 2 



(35) 



In addition, the metric functions of the B region decay 
exponentially as the horizon is approached, as demon- 
strated in Sect. 11111 That is, 



where k_ is the surface gravity of the Cauchy horizon; 
a positive quantity that is defined in Eq. ( p37| (with m 
replaced by m B ). The negative sign in front of the ex- 
ponent indicates that we are approaching the Cauchy 
horizon from between the horizons, and so X B will be 
negative as this approach takes place. This fact is not 
as obvious as it may seem for the loop black hole. The 
only reason it is correct is because ms is fixed; this sub- 
tlety will be discussed later in this section and in the 
Appendix. 

In the same way that mc~m B will be the energy of the 
ingoing shell, the mass difference mo — ms will be equal 
to (the negative of) the energy of the outgoing shell, at 
least in a classical setting. That is, the outgoing shell has 
negative energy density and we have m B > mg. How- 
ever, this condition results from the Einstein field equa- 
tions, and in the case of the outgoing shell we cannot 
necessarily make the same conclusion because we don't 
have the extension to the outside of the black hole that 
we do for the ingoing shell (i.e. the emission of the outgo- 
ing shell from the collapsing star occurs inside the black 
hole). This means that without knowledge of the effective 
field equation we can not know with certainty whether 
mo will be greater than or less than m B . However the 
loop black hole admits a maximally extended universe, 
suggesting that m B > m B in the loop black hole, as in 
the classical case. 

Fortunately, whatever the case, all we need to know 
is that to B ^= m B , because this is what produces the 
separation between the apparent and Cauchy horizons 
that is necessary for mass inflation. We will soon see 
explicitly why we don't need to specify the direction of 
mass change across the shell. 

From Eq. ( 32 ) we have 



\*a\ = 



Xc 
X B 



\Xj 



(37) 



Using this we are going to examine what happens to Xa 
when we take v — > oo, and we will find in fact that it 
grows unboundedly large. This is not obvious because 
while X B approaches zero as v — ¥ oo, so does Xc (after 
all we have X B — Xc at v = oo). The point is that 
X B decays faster than Xc as the Cauchy horizon is ap- 
proached, and so the limit of Xc/X B will diverge. More 
explicitly, given mc = m B + v 1 we have that 



X c 



X B oc 



(36) 



(r - 2m c )(r - 2m c P 2 ) 

(r + 2m c P) 2 
(r - 2m B ){r - 2m B P 2 ) - 2y(r)v~"' 
(r + 2m B P)(r + 2m B P + APv^) 
(r + 2m B PfX B - 2y(r)y-"i 
(r + 2m B P) 2 ' 



for large v. Here we have used 

y(r) = (1 + P 2 )r ~ Am B P 2 



■0(v- 2 ~<) 
(38) 

(39) 
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Thus, as v — > oo we see 
Xc 



■»(r-)»" 1 V 



y(r_)^e K ^, 



(40) 
(41) 



which indeed diverges. 

The assumption that the outgoing shell produces a 
mass difference Tod / means that is nonzero 
on the Cauchy horizon, and so by Eq. (37) we conclude 
that Xa diverges as 



\X, 



(42) 



As stated above, this is independent of whether mjj is 
greater than or less than vhb since this only acts to 
change the sign of Xd on the Cauchy horizon. Thus, we 
see that indeed there is divergence at the Cauchy horizon. 
In addition, this is the exact same form of divergence that 
is observed in classical mass inflation. Note that this di- 
vergence is also similar to the form that was observed 
from the perturbation calculation in Sect. (Ill I. 

As a quick aside, we find that y(r_) = j/(2tobP 2 ) = 
4m B P 2 (P 2 - 1) < since P < 1. This tells us that 
Xc/Xb < for large v which, since Xb remains nega- 
tive, tells us that X c is positive as the Cauchy horizon is 
approached. Recall our earlier remark noting that if the 
loop black hole is subjected to ingoing radiation the met- 
ric function need not remain negative for large v before 
the Cauchy horizon is reached. This is a result of the fact 
that, unlike the Reissner- Nordstrom black hole, the in- 
ner horizon of the loop black hole inflates with increasing 
mass. This produces a strange effect where, when con- 
tinuous ingoing radiation is present, an outgoing geodesic 
can pass by the radius at which the Cauchy horizon will 
be located. The geodesic will then reach the apparent 
horizon, turn around, and straddle the apparent hori- 
zon before settling on to Cauchy horizon at v — oo. We 
discuss this situation more thoroughly in the Appendix. 
The positivity of Xc observed here is the shell model 
analog of this phenomenon. This has no consequence on 
the divergent behaviour of Xa- 

We have seen that Xa, and therefore the metric 
function Fa diverges at the Cauchy horizon. Similar 
reasoning yields the same conclusion for the Reissner- 
Nordstrom black hole, except that it is the function 
Ja = 1 — 2m a/ r + e 2 /r 2 that diverges, thus indicating 
that the mass rriA diverges, hence the name mass infla- 
tion. In the loop black hole, however, Fa — > oo cannot 
be achieved from a diverging mass parameter (see Eq. 
(31)). Indeed, Fa remains perfectly finite for uia — > oo, 



and so referring to our result as mass inflation is arguably 
a little disingenuous. However for historical reasons and 
because our result follows from the same physics and has 
the same characteristics as classical mass inflation, we 
will continue to call it as such. 

It is not hard to see that the only way for Fa to actually 
diverge is for the mass to take a negative value of uia —> 



— r_/2P. It is not clear how to interpret this strange 
result. We interpret the parameter m in the loop metric 
as being the mass because at infinity Madm = (1 + 
P) 2 m, where Madm is the classical ADM mass. To the 
future of the outgoing shell, however, the parameter m 
no longer has this classical interpretation. Alternatively, 
one could argue that we should define m in the context 
of a quasi-local mass computed from F. However this 
interpretation is based on the Einstein field equations 
|26j . and thus the same interpretation can certainly not 
be applied inside the loop black hole. We therefore simply 
consider m to be a parameter of the spacetime, and so 
are not concerned with the possibility of negative tua, 
since this negativity occurs in a region of spacetime where 
there is no good reason to interpret uia as being a mass 
anyway; the result is simply a geometric requirement. 

When the limit tua —> — r_/2P is taken at the Cauchy 
horizon the Ricci and Kretschmann scalars computed 
from the loop metric go as R ~ (r_ + 2m aP)~ A and 
K ~ (r_ + 2rriAP)~ 6 respectively, thus indicating that 
they diverge in this limit. That is, the phenomenon seen 
to occur here causes a scalar curvature singularity at the 
Cauchy horizon, the same result as that obtained in clas- 
sical mass inflation. However the DTR analysis we have 
performed only tells us what happens to Fa- Since for 
the loop metric Fa is a quadratic in uia there are two 
possible solutions for tua, yielding an ambiguity in how 
rriA approaches this negative value as, say, a function of 
v. In the Reissner-Nordstrom black hole there is no am- 
biguity since the metric function is linear in the mass. 
In our case all we know is that Fa increases, eventually 
diverging at the Cauchy horizon. Solving for mA(FA) we 
find that as Fa increases one of the solutions goes nega- 
tive and approaches — r_/2P continuously. In the other 
solution uia first goes to positive infinity and then jumps 
to negative infinity before approaching ~r_/2P from the 
negative side. 

To understand how these two solutions arise we plot 
at fixed r in Figj4] Fa as a function of mj, where P 
and ao are some fixed values. We see that, given some 
initial value of Fa (presumably negative) , there are two 
possible values of m^; either on the left or the right side 
of the minimum in Fa- The DTR relation tells us that 
Fa increases as the Cauchy horizon is approached, and 
we see from the figure that as Fa is increased the value of 
ttia will either decrease or increase depending on which 
option we choose for its initial value. If it is to the left of 
the minimum in Fa its value will decrease continuously 
and eventually reach —r/2P as Fa diverges. If it is to 
the right of the minimum, however, its value will increase 
and diverge to infinity when Fa — r 4 /(r 4 + ctg) before 
reemerging at negative infinity and approaching —r/2P 
from the negative direction. 

Note that, of course, Fa is also a function of r and 
so the plot will change as the evolution occurs since the 
value of r will also be changing. We find that the same 
qualitative behaviour occurs in the full scenario as well. 
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That is, along an outgoing null geodesi c r g oes like r(v) 



r_ + k_ e K - v for large v 
be substituted into Eq, 



(38) 



see Sect. 

When" 
ved, we find that the two 



III), and this can 
this is done and 



the resulting quadratic is so 
solutions follow the behaviour described above. 

Thus, we see that which evolution of vtia actually takes 
place depends critically on the initial condition placed 
on rriA- Unfortunately within the confines of our model 
we have no way of specifying what this initial condition 
should be. This is because the power law decay given by 
Eq. (34) is only valid in the large v limit. Presumably 



there will be a nonzero influx of radiation during early 
times which we have no knowledge of, and so vtla will 
have already gone through some evolution by the time 
the large v regime is reached. Despite this limitation, 
however, we can perform a simple thought experiment in 
which we suppose that there is some point in the past, 
shortly after the collapse of the star, at which the influx 
starts. That is, we imagine that for some v = vq the 
energy density of the ingoing radiation transitions con- 
tinuously from zero to nonzero. In the thin-shell picture, 
this initial ingoing radiation will intersect the outgoing 
radiation shell at some r = tq, producing an intersec- 
tion two-sphere located at So — (vo,ro). In this scenario 
we have, for example, that mc(v = Vo) = "niB because 
the black hole has only just begun to accrete extra mass. 
In addition, it is clear from physical considerations and 
the assumption that wia should evolve continuously that 
we also have m,A{v — vq) = mp, thus giving us an initial 
condition for vtia- From here, in order to determine which 
evolution of wia will ensue we note that the minimum 
in Fa seen in Fig. [4]) is located at mass vtia = r/2P. 
Therefore, we see that if m^^o) = mjj > r$/2P then 
wia will follow the evolution that takes it to infinity and 
back around. Alternatively, if mr> < tq/2P then itia 
will proceed to its final value of — r_/2P continuously. 
Given the scenario in our thought experiment, the evo- 
lution of ttla is thus entirely dependent on the value of 
m b ) whether mo is greater than or less than mj and by 
how much, and the radius r$ at which the cross streams 
first intersect. 

Note that for m — > ±oo the metric functions approach 



F -> r 4 /(r 4 + o§) and G -> 16m 4 P 4 /(r 4 + a§) -> +oo, 
yielding a diverging volume element \J—g. However both 
the Ricci and Kretschmann scalars, R and K, remain 
perfectly finite for diverging m, a fact that remans true 
when we include a functional dependence m = m(v); i.e. 
the diverging derivatives of m(v) do not spoil the fact 
that R and K are finite. 

Furthermore we find that there are no nonscalar sin- 
gularities in this limit. The presence of a nonscalar sin- 
gularity is signalled by the divergence of at least one of 
the tetrad components of the Riemann curvature ten- 
sor when the tensor is evaluated in a parallely propa- 
gated orthonormal frame (PPON frame) |27j . If all of 
the components remain finite this means that the ob- 
server traveling along the geodesic defining the PPON 



frame will not experience any diverging tidal or inertial 
forces upon approaching the point in question. Given a 
radial geodesic U a , we construct a PPON tetrad by let- 
ting the first vector be tangent to this geodesic e" = U a 
and, given spherical symmetry, we allow two of the other 
vectors e% and to lie in the two-sphere, at which point 
ei follows directly from orthonormality r/ij — e°^e^g a p. A 
radial geodesic (inside the black hole) on the loop metric 
takes the form 



E 
G' 



jjr = _ 



G 



(E 2 -G), f/ e = [/* = 0, (43) 



where E is a constant (see Sect. III). Thus, in {t,r,9 



coordinates one finds the PPON frame generated by this 
geodesic to be 




-(£2-G),0,0 



(44) 



0,0, 



.0 



0,0,0, 



1 



H sin 9 



Projecting the Riemann curvature tensor on to this basis, 
Ra/3~f8 we find, somewhat miraculously, 



R 



ijki 



c,' oi e 



that all components remain finite when m is taken to 
infinity. 

These results suggest that taking the limits m — > ±oo 
does not result in a curvature singularity. Note that, 
like the scalar invariants, the curvature appears to re- 
main finite even when we include a functional dependence 
m = m(v) (note that the tetrad was evaluated with re- 
spect to (v, r, 9, tfi) coordinates when performing this cal- 
culation). Of course in this dynamical case the tetrad in 
Eq. ( |44[ ) is no longer a PPON frame, and unfortunately 
finding a closed form for radial geodesies in this space- 
time is nontrivial since t is no longer a Killing direction. 
Even so, when m = m(v) then Eq. (44) still represents 
the frame of an accelerated observer, and the fact that 
the tetrad components Rijke remain finite when m(v) 
and its derivatives diverge means that, at least for this 
accelerated observer, the spacetime curvature appears to 
remain finite. Note that Rijke also remain finite if we 
use, for example, the static orthonormal tetrad given by 
ef = (1/VG, 0,0,0), e£ = (0,^,0,0), along with 
and e J given as they are in Eq. ( 44 ) . 



To summarize, we have been unable to find a way in 
which the loop metric becomes singular when m — > ±oo, 
apart from a divergence in the volume element, a coordi- 
nate artifact. This, coupled with the fact that the metric 
functions F and G take the same value when m is both 
positive and negative infinity (since G cx m 4 for large ±m 
, then G attains the same limit for either sign), suggests 
that evolution of tua through infinity to — r_/2P may in- 
deed represent a perfectly well behaved and continuous 
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evolution of spacetime. Note that ttia is a parameter of 
the spacetime that is only accessible deep inside the black 
hole and is not observable from the outside. Further- 
more, we can't identify mj with any sort of quasi-local 
mass since the Einstein field equations are not expected 
to hold in this region of spacetime. We thus identify ttia 
only as being a parameter of the internal spacetime, one 
whose evolution yields a continuous sequence of space- 
times even though the observed evolution of rriA itself is 
discontinuous. 

Whichever solution ttia follows, it will always limit 
to the value of — r_/2P as the Cauchy horizon is ap- 
proached, and this certainly is a surface of infinite cur- 
vature. It is interesting to note, given that the spacetime 
appears regular for m — > ±00, that the DTR analysis 
ferrets out the only way that the loop metric knows how 
to become singular. 

The model used and the DTR analysis performed here 
are quite general. Any black hole with a Cauchy horizon 
- classical or quantum - will be hard pressed to avoid the 
same result. 

The two solutions observed for here for the evolution of 
m are also seen to occur when this problem is analyzed 
using the Ori model. We discuss the implications and 
limitation of this approach in the Appendix. 



VI. CONCLUSIONS 

We have studied the question of whether or not the 
phenomenon of mass inflation occurs in the loop black 
hole. In attempting to provide an answer we employed a 
simplified model in which the ingoing and outgoing null 
radiation (whose counter-streaming classically results in 
mass inflation) are regarded as null thin-shells. Within 
this framework, the generalized DTR relation allows us 
to make conclusions about how the metrics in the four 
resulting spacetime quadrants match at the intersection 
two-sphere. Using the DTR relation is of particular util- 
ity in our case because said matching does not directly 
rely on Einstein's field equation; this is good because we 
are attempting to work within a loop quantum gravita- 
tional framework as much as possible. 

There are four principle assumptions that were made 
during our analysis: 

i) The LBH metric is a valid effective metric de- 
scribing a black hole spacetime. 

ii) The spacetime is no worse behaved on the two- 
sphere S than it is across the shells. Equivalently, for 
every point on S there is a local coordinate system in 
which the metric is continuous and piecewise continu- 
ously differentiable. 

iii) Each spacetime quadrant bounded by the colliding 
shells is itself described by an effective loop metric of the 
form ([!]). 

iv) The energy of the ingoing shell as a function of 



v follows the same inverse power law that is known for 
continuous ingoing radiation. 

The conclusion of our analysis is that the Cauchy hori- 
zon for the LBH remains unstable, implying that the loop 
black hole is not a fully adequate effective metric for de- 
scribing expected black hole corrections by loop quantum 
gravity. The DTR relations force the curvature to diverge 
as the parameter tua — > — r_/2P, a fixed negative value 
dependent on the parameter of the exterior solution. We 
find that rriA can approach this value either directly or by 
passing through infinity; in the latter situation we have 
been unable to uncover any sort of curvature singularity 
as rciA — > ±00. 

We demonstrate in the appendix that our results are in 
part corroborated by the Ori model |16j . Apart from rely- 
ing directly on Einstein's equations, this model is limited 
insofar as we are unable to specify an equation of state 
between the pressure and energy density of the outgo- 
ing shell for the LBH. However in the limit the pressure 
is negligible the model predicts the same behaviour that 
was discovered using the DTR model. 

The analysis performed here will be applicable to any 
other quantum gravitationally corrected black holes that 
yield effective spacetime metrics with inner horizons, 
such as found in PH5]. Indeed model calculations of the 
type performed in section II for non-commutative black 
holes [28] indicate that their interiors are even less stable 
than that of the LBH For them to avoid the same 
fate as the LBH will require very special properties of 
the effective metric. Either they will not contain inner 
horizons due to their own 'quantum repulsion' (a diffi- 
cult feat which is seemingly impossible for black holes 
with angular momenta) or these inner horizons will have 
some kind of cutoff property that allows metric functions 
in the B region to avoid the behaviour in ( 36 ). It will cer- 



tainly not be adequate to placate the r — singularity, 
the primary goal in deriving these quantum black hole 
solutions, to attain a purely regular spacetime. 

Finally, we wish to ponder whether or not the dis- 
covery of an unstable Cauchy horizon should be an en- 
couraging one or not. The primary goal of developing 
quantum gravity black hole solutions seems to be to reg- 
ulate the singularities present in their classical counter- 
parts. Of course the structure of spacetime at the thin 
shell juncture and the Cauchy horizon will undoubtedly 
be modified by quantum gravitational effects, whatever 
the correct theory of quantum gravity is. Stabilizing the 
Cauchy horizon should therefore be regarded as one of the 
major goals of quantum gravity. However a traversable 
Cauchy horizon yields a spacetime (and a theory) that 
is no longer deterministic. Our results indicate stabiliz- 
ing the Cauchy horizon remains a major challenge for all 
quantum gravitational theories. 
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Appendix: The Ori Model 



For the LBH the situation is more complicated. To see 
this, consider the metric 



ds 2 = -G(r, v)dv 2 + 2* ffi V \ drdv + H(r)dn 2 , (47) 



where 



G(r,v) = 



F{r,v) = 



F(r, v) 



(r - 2m(v))(r ~ 2m{v)P 2 )(r + 2m(v)Pf 
(r - 2m{v)){r - 2m{v)P 2 )r 4 



(r + 2m(v)P) 2 {r 4 + a 2 ) 



H(r) = r 2 + 



2 , "0 



(48) 



In this appendix we aim to outline how the Ori model 
of mass inflation [IB] applies to the loop black hole and 
why it fails to be fully predictive in this case. Even so, 
we will see that with a simple assumption it nevertheless 
gives the same result as was reached by using the DTR 
model. 

In the Ori model outgoing radiation is modelled as a 
null thin-shell, but the ingoing radiation is continuous. 
In this sense it is intermediate between the full model 
(in which both ingoing and outgoing streams are contin- 
uous), and the DTR model that considers both streams 
as null thin-shells. We will start by considering the loop 
black hole in the presence of continuous ingoing radia- 
tion and how this changes the spacetime structure. Un- 
like the analysis in the main body of this paper, we will 
be using Einstein's equations to make conclusions about 
curvature. 



A. The Vaidya loop metric 

Recall that the advanced time v is defined as v = t+r*, 
where r* satisfies 



dr* _ 1 



(45) 



and the metric functions G(r) and F(r) are defined as 
in Sect. [TTJ If we cast the loop metric into (v, r, 8, <f) 
coordinates it takes the form 



ds 2 = -G(r)dv 2 + 2a 



I G(r) 
F(r) 



dvdr + H(r)dn 2 , (46) 



where the warped function is H(r) = r 2 + a^/r 2 . Null 
radiation entering a Reissner-Nordstrom (RN) black hole 
as a continuous stream takes the form of null dust with 
stress-energy T a p = fx(r,v)(d a v)(dpv). The net effect of 
this is to change the mass parameter from a constant to a 
function of v. m — > m(v), with the energy density taking 
the form /x(r, v) = m! (v)/4irr 2 . 



which we will refer to as the Vaidya loop black hole. We 
see that the inner and outer apparent horizons occur at 
r = 2m(v)P 2 and r = 2m(v) respectively. The Einstein 
equations then imply 



ydyn — 



(1 + P) 2 r 2 (r 4 - al)(r - 2m(v)P) 
47r(r 4 + a 2 ) 2 (r + 2m{v)P) 



mf(v), (49) 



T ™ = 2,(r + 2m(v)PY m {v) ' T » = Sm ^ ' 

where the "dyn" stands for dynamical; we include it to 
emphasize that this is the stress-energy in addition to 
that which is already present for the original "static" 
LBH. We will call the stress-energy of the original metric, 
which is nonzero, TS^ at ; its exact form is unimportant. 



Clearly, the stress energy in Eq. (49) does not cor- 
respond to null dust. For oq = P = it reduces to 
the usual form of T^ n = m! (v) / Anr 2 . In addition, for 
large r (where we know that the black hole is nearly 
Schwarzschild and classical theory should approximately 
hold) it reduces to T^ n = (1 + Pfm! (v)/4Trr 2 , T^ D = 
—Pm'(v)/2-K which, since the ADM mass Madm 

1 + P) 2 m, shows us that the energy component T^ n 
reduces to the usual form that we expect and that Tg@ n 
is of order P, which is assumed small. 

Ultimately, the strange form of the stress-energy re- 
sults from the failure of the Einstein equation to remain a 
valid field equation deep inside the black hole. Note that 
the ingoing energy flux, which is proportional to T^ n , 
becomes negative in between r = 2m(v)P and r = y/oQ. 
If we interpret the modified stress-energy as being an ef- 
fective quantum gravitational correction, then this neg- 
ativity represents the correction that allows avoidance 
of the central singularity. Quantum gravitational effects 
work against the incoming radiation such that near the 
center of the black hole the flux becomes negative, thus 
avoiding collapse to a singularity. As the radiation con- 
tinues to smaller r it enters the second asymptotically 
flat region and the flux becomes positive again, where 
it is then emerges as radiation from the corresponding 



14 



white hole. Regardless of the interpretation, the fo- 
cusing theorem guarantees that the null energy condi- 
tion must not hold for null radiation that undergoes the 
"bounce" evolution characteristic of the warped factor 
H(r) = r 2 + a^/r 2 . Hence the form of T^ n is to be 
expected given the nature of the loop black hole. 

A valid criticism of the preceding discussion is that, 
since the loop black hole was constructed as a model for 
the quantum gravitationally corrected vacuum solution, 
we should not presume to simply replace the parameter 
m with m{v) in the presence of ingoing null radiation. In 
response, in the 2 nd paragraph of Sect. [v]we argued that 
m should be constant within a given spacetime quadrant. 
This reasoning is independent of how many ingoing null 
shells are present in the spacetime. Conceptualizing the 
continuous radiation as a sequence of "strips" , we expect 
m to be constant within a given strip. The Vaidya loop 
metric ( 47 ) is the limit in which we have infinitely many 
ingoing null shells, each of infinitesimal energy. It fol- 
lows that, given a constant m within each infinitesimally 
thin strip between adjacent shells, this is equivalent to m 
being constant along lines of constant v. 

In the large v limit we expect the energy density of the 
ingoing radiation to fall off as an inverse power, as shown 
by Price [15] ■ This results from backscattered gravita- 
tional radiation, the origin of which is the surface of the 
collapsing star. Since this process occurs outside of the 
black hole, where classical theory is expected to hold to 
good approximation, we can be confident that the same 
result will apply in our case as well. This means that for 
large v we have 



m(v) = M — large v, 



(50) 



where M is the final mass of the black hole and we ex- 
pect fj,(v) oc u 1 ^ 7 , where 7 > 12. Note that the Cauchy 
horizon is located at coordinates r = r_ = IMP 2 and 
v = 00. 

In addition to continuous ingoing radiation there is 
also a null thin-shell that acts to separate the spacetime 
into two regions, one to the future and one to past of the 
shell. Consider the behaviour of outgoing null geodesies 
in the black hole, specifically in the presence of ingoing 
radiation. To this end we define a function 



f(r,v) = y/G(r,v)F(r,v) 

(r — 2m(v))(r — 2m{v)P 2 )r 2 
r^Ta 2 • 



(51) 



We see from Eq. (471 that the two radial null directions 



are given by setting dv — for the ingoing direction and 
2dr = fdv for the outgoing one. Hence 



dr = / 
dv 2 



(52) 



is the equation obeyed by the geodesic r(v) for an out- 
going photon. 



Inserting m(v) — M — fi(v) in to Eq. (p3lj) yields 



where 



f(r, v) = f(r, 00) + g(r)fi(v), 



/(r,oo) 



9{r) = 



(r-2M)(r-2MP 2 )j 



2[(1 + P 2 )r- 4MP 2 }r 2 



(53) 

(54) 
(55) 



and we have neglected a term proportional to ^(v) 2 (this 
is justified since we are working in the v — ¥ 00 limit). Ex- 
panding f(r,v) around the Cauchy horizon r_ = 2MP 
(which is the apparent horizon in the v — > 00 limit) gives 

f(r ~ r-,v) ~/(r_,oo) + /' (r_,oo)(r - r_) 

+ g{ r -)^(v) + g'(r-)(r - r-)(i(v). (56) 

The first term is zero (by definition) and we neglect the 
last term because it is 2nd order (since r — r_ is also 
vanishingly small). We identify 



;/'(»•-, 00) 



4M 3 P 4 (1 - P 2 ) 
16M 4 P 8 + ak 



(57) 



as the surface gravity at the Cauchy horizon. Note that 
the surface gravity in the LBH is not in general given by 

\\d r ^/GF\ but rather by k = | \d r G\. However at the 

horizons these two expressions happen to give equivalent 
answers. Evaluating g(r ) in terms of r_ = 2MP 2 gives 



2{P 2 

9{r-) = 4T 

r_ 



l)r 3 



(58) 



which is a negative quantity. 

Putting this together with the null geodesic equation 
dr/dv = f/2 yields an ODE for r(v) in the large v limit 



dr 



dv 

whose solution is 

K_(r — r_) 



+ K-(r — r_) 



(P 2 - l)r 3 _ 



(59) 



(60) 



{P 2 



1 v 3 

— E 

71=0 



ft, 







d n ^{v) 
dv n 



where a is a constant. The summation term represents 
the correction to K_(r — r_) from the v dependence; in 
the static case we only have the exponential term. In 
the v — > 00 limit we expect the exponential term to be 
negligible compared to the sum, and therefore will ignore 
it from this point on. Putting this result back into the 
expansion of / gives 



f{r 



(P 2 - l)r 3 



E 



1 



d n n{v) 
dv n 



(61) 
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This is the value of / that an outgoing photon will ob- 
serve for large v as it travels along its geodesic. 

Here we encounter an apparent mathematical contra- 
diction. We expect that the LHS of Eq. ( 60 1 should be 
positive since k_ is defined to be positive and r > r_ 
before the geodesic reaches the Cauchy horizon - yet we 
have shown here that the RHS is negative! (For the RN 
black hole the result is exactly the same except that the 
factor in front of the sum is l/r_ , which is positive). 
Similarly, we see that / is positive, which seems to sug- 
gest that for large v the outgoing photon is no longer 
between the horizons - by Eq. ( 52 ) , that its radial value 



r is increasing! 

The problem arises in our assumption that r > r_ for 
large v. The key difference between the RN and LBH 
cases is that in the presence of ingoing matter the appar- 
ent inner horizon of the RN black hole deflates, meaning 
that an outgoing photon has to "chase down" the inner 
horizon before finally intersecting it at v = oo, where the 
inner apparent horizon coincides with the Cauchy hori- 
zon. In the loop black hole, however, the inner apparent 
horizon inflates for increasing v, as is obvious by its form 
r = 2m(v)P 2 . This means that, rather than receding, 
the horizon instead rushes up to the meet the outgoing 
photon. This results in the photon first passing by the 
radial value that the Cauchy horizon will be at (resulting 
in r < r_ ) before penetrating the inflating apparent hori- 
zon (resulting in / > 0), turning around and straddling 
the horizon before again reaching the horizon at v = oo, 
which becomes the Cauchy horizon. 

This behaviour can be seen in Fig. [5]) where we have 
numerically solved r(v) for general v (i.e. we use the full 
form of Eq. (51) in Eq. ( 52 ) ) but for simplicity contin- 



ued to assume an inverse power law for fj.(v). Specifically 
we used = (k-v)^ 1 ; this choice is fairly unphysical 
since the inverse power should be much higher, but this 
lets us see the behaviour much more easily. The same 
qualitative behaviour occurs for more realistic fall offs as 
well. The parameters used in the solution are M = 1, 
P 2 = 0.1 and Oq = 0.15, giving the radius of the Cauchy 
horizon as r_ = 0.2. All of the stages in the evolution de- 
scribed above are clearly visible in the plot. Despite this 
strange and unexpected behaviour, the geodesic evolu- 
tion remains perfectly well defined and the Cauchy hori- 
zon is still a surface of infinite blue-shift. 



B. The null thin-shell formalism 

We review here the null thin-shell formalism [231 [H] , 
outlining the results relevant to the situation at hand. We 
note that the results of this section rely directly on Ein- 
stein's equation, and so to use it within a loop quantum 
gravitational framework is at best a first approximation. 

We wish to determine the relationship between two 
metrics g^ with corresponding coordinates a;™ on either 



side of a thin null shell S. Suppose that E is equipped 
with some coordinate set (A, A ), where A is a parameter 
which runs along the null generators of the shell and A 
runs over the remaining two coordinates. The surface 
E is described by x±(y a ), with tangent vectors e± a — 
dx±/dy a ; from these we identify the vector tangent to 
the generators k a and those transverse to them e A 



k a 



dx a 



dx a 
dQ 1 ' 



(62) 



For simplicity we have removed the ± notation; this 
should not introduce any ambiguity, keeping in mind that 
k is a vector given by k — d/dX, and the only thing that 
changes between the two sides of the shell is the coordi- 
nates in which we evaluate this vector. 
Since this is a null shell, we have 



k^k — — k(%e A 
and the induced metric is given by 

_ a P 

o~ab = 9 a pe A e B . 



(63) 



(64) 



In order that the shell itself have a well-defined intrinsic 
geometry the induced metric should be independent of 
which side of the shell it is being evaluated on, implying 



[g-ab] = 0. 



(65) 



To complete the basis over all spacetime dimensions 
we need to include another vector that isolates the part 
of the metric g a p transverse to the shell. To this end, we 
introduce an auxiliary vector N a that solves 



N a N a = = N a e% 



N n k a = -1. 



(66) 



These conditions specify N a completely. Note, however, 
that one needs to compute this vector separately for each 
side of the shell. With this, the total metric can be shown 
to be given by the completeness relation 



g~- = -k a N? - N a k fj + a AB e%e B , (67) 

where a AB is the matrix inverse of cab (note that this 
is just a special case of Eq. (27)). 



We also introduce what is known as the transverse cur- 
vature Cab = \i.^-Ngap)e-Z e b i w hich more explicitly is 



Cab 



-N a e 



(68) 



where Latin indices run over three coordinates: y a 
(X,9 A ). The stress -energy of the shell is |24j 



Tf = (-k^r'S^Sir), 



(69) 



where is tangent to an arbitrary congruence of time- 
like geodesies parameterized by r in which the geodesies 
cross E at r = 0. The 8{t) factor represents the assump- 
tion that the shell is infinitely thin, and is equivalent to 
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the divergence of the Einstein tensor mentioned above. 
The tensor S a ^ is given by 

S aP = iik a k p + j A (k a e A + elk 13 ) + pa AB e\e p B , (70) 

where the energy density, current and surface pressure 
are given by 



1 



AB 



Cal 



3 



1 

8^< 



AB 



[C 



AB 



AA I • 



(71) 



Note that the factor (— fc^u^) -1 means that T^f is ob- 
server dependent; a measurement of the shell's stress- 
energy will depend on the geodesic followed by the ob- 
server doing the measuring. We can thus define the 
"physical" quantities measured by this observer: 



Mphys — ( k^U^) 1 fl, jphy S — ( k^U^) 1 j A , 



Pphys = {-k^) *p 



(72) 



AA 



N a k°!gk" = k, where n is the " 



ac- 



Noting that C. 

celeration" of the shell, fc Q fl fc' 9 = Kk a , we obtain from Eq. 



( 71 ) the matching condition 

[«] = — 87175, 



(73) 



which, via Raychaudhuri's equation (and making use of 
hypersurface orthogonality), is equivalent to 



[k]6 = 87r[T Q ^fc Q fc /3 ], 



(74) 



with 9 the expansion of S (with respect to the generating 
parameter A). T a p is the stress-energy of the surrounding 
spacetime, not of the shell. 

We close this subsection by considering what happens 
when we choose a different generating parameter, A, to 
describe the shell. Under the reparametrization 



A — > A(A, 9 ), 



wc define 



9 



implying that the generating vector transforms as 



k a 



<J 



l k a 



(75) 



(76) 



(77) 



and, less obviously, the stress-energy quantities defined 
in Eq. ( 71 I transform as 



fi = gjj, + 2c A j A + c A c A g 1 p, 
J A =3 A + c A g- 1 p, 
P = .9~V 



(78) 



In total it can be shown that the stress-energy transforms 
as S a/3 = g~ 1 S al3 and therefore that the combination 
{—k ll u^)~ 1 S a P is invariant under the the reparameteri- 
zation. That is, we have invariance of the total stress- 
energy of the shell T^f . 



C. The evolution of m+(v) 

Taking ingoing radiation to be continuous, giving us 
the Vaidya spacetime Eq. (47), and outgoing radiation 
modeled as a null thin-shell S that acts to split the space- 
time into a future region (+) and a past region (-), gives 
a scenario displayed in Fig. [6j 

Since on either side of the shell only an influx of radi- 
ation is present, it follows that the spacetime should be 
described by a Vaidya loop metric on both sides, given 
by 



ds\ = —G±dv± 




drdvA 



Hdfl 2 



(79) 



where the metric functions G± and F± are defined as 
above, but with m(v) replaced by m±(v±). Henceforth 
for convenience we will drop the subscript on the ad- 
vanced time on the past side w_ = v since we expect 
that here the spacetime behaves as outlined in part A of 
this appendix. The mass will be given by the function 
m_ (v) = M — fi(v) and the coordinate v will still deter- 
mine the location of the Cauchy horizon (at r = 2MP 2 
and v — oo). 

To investigate the behaviour of m + , we begin by not- 
ing that the tangent vectors to the shell transverse to the 
generators are 



dx a 



(0,0,1,0), 



The induced metric is thus 



o~ab 



dx a 



1 



H ( r ) In -2 
w \0 sm 



(0,0,0,1). (80) 



(81) 



from which the condition specified by Eq. (65), [(tab] = 
0, gives us [H(r)] = 0. The solution that we take 
from this condition is that the coordinate r is continuous 
across the shell. 

With the knowledge that r is continuous across E we 
can immediately read off the geodesic of the shell as seen 
from either the past or future side (i.e. in either coor- 
dinate system), and in so doing also obtain a relation 
between v and v + (which is valid only along E) 



2dr = f + dv+ = f-dv, 



along E, 



(82) 



where for convenience we have defined 



f± = VG^F± = ( r - 2m ±( v ))( r - 2 2 m ±( v ) p2 y (83) 



It it important to keep in mind that in all subsequent 
equations the coordinate r is actually a function of v 
since everything is evaluated on the shell. Specifically, 
the form of r — r(v) is given by Eq. (60) and looks like 
the solutions plotted in Fig. [5] 
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Defining the tangent vector k a to be 



k a = 



dx a 
dv 



we find 



dx a _ 
dv 



f dx a 

(1,^,0,0), K = ^± 



(84) 



on the past and future side of shell, where — 
(v±,r,9,(f>) and v_ = u. We wish to compute the accel- 
eration k — —N a k a gkP on either side so that we can use 
the equation [k] = —8irp. In order to do this we need to 
first compute the auxiliary vector N a ; the form of which 



is fully determined by the conditions in Eq. ( 66 1 . Unlike 



k, N a differs on the two different sides of the shell. On 
the past side the conditions on N a yield 



JVf! = (0, 




,0,0) 



and on the future side they give us 



A^ = (0,-^,0,0). 



(86) 



(87) 



Computing k = —N a k a pk@ on both sides of £ and 
noting that m' ± (v + ) = j^m' ± (v) yields the differential 
equation 

r 3 (l + P) 2 (r -2m+P) 1 dm+ 

(r 4 + al){r + 2m + P) j+~ dv 
r 3 (l + P) 2 {r - 2m_P) 1 dm_ 

(r i + a 2 ){r + 2m_P) JZ~dV 



P 



(r + 2m_P)(r + 2m+P) 



/_ (m_ — fn+) 



where m_ = m_ (v) = M — fi(v) is a known function and 
r = r{v) is the solution of r'{v) = f-/2 (as per Eq. (60 1). 



From this the metric functions are purely functions of v: 
f± = f±(r(v),m±(v)). This is because we are evaluating 
everything on the shell, which follows a known geodesic. 

Note that we could just as easily have done this anal- 
ysis using r to parameterize £ instead of v. Replacing 
v — > r then p — > (dv/dr)p = (2/f-)p which, when sub- 
stituted into Eq. (88 yields the same equation but with 
the RHS multiplied by /_/2. 

The only thing left to do before examining the be- 
haviour of 77i+ is to specify the pressure p. Unlike the 



RN case we cannot assume that the shell is pressureless 
due to Eq. (49). We therefore need some other equa- 
tion of state (EoS) p — p(p), where p is the energy- 
density of the shell, to specify the pressure. Using Eq. 
( 71 ) we can compute the energy density of the shell 
^ :<j AB [Cab\, yielding us some function of r, m_ 
which, given p = p(p), can be inserted into the 



M = "Si 
and T7i + 

RHS of Eq. (|88| to give us an equation without any 



(g5^ unknown, free parameters. 



There are two equivalent ways in which we could try 
to obtain an EoS given our model. The first is to look at 
the stress-energy of the ingoing radiation, Eq. (49), and 
attempt to simply read off the appropriate EoS. As can 
quickly be seen, however, this simply is not possible given 
the form of the stress-energy we have; an EoS does not 
exist. Alternatively, we could compute the energy density 
H and pressure p of the shell E (using Eq. (71)) in the 



case that there is no ingoing radiation. For example, if we 
compute the pressure in this scenario we obtain Eq. ( 88 ) 
but with the first two terms set to zero. However this 
entails the same problems as the first approach, because 
the scenario of the outgoing shell in an otherwise static 
black hole is really just a special case of an outgoing 
Vaidya spacetime, and so the stress-energy of these two 
situations will be of the same form. 

Without being able to specify the pressure of the shell, 
the Ori model for the LBH loses its predictive power; 
each choice for the p will give a different evolution for 
to_|_(i>). Despite this, it is interesting to note that for the 



choice p = the solution m + (v) of Eq. (88) displays the 
same behaviour as that obtained for itia(v) in the DTR 
analysis. That is, we find that if the initial condition 
for m_|_ is small enough then it evolves continuously to 
a value of — r_/2P, and if the initial condition is higher 
then the solution is seen to diverge to infinity for finite 
v. These are exactly the two options found for m^. 



To see this explicitly, we numerically solved Eq. (88) 
in the large v limit. We set the parameters to M = 1, 
P 2 = 0.25 and a§ = 0.1; this gives -r_/2P = -MP = 
—0.5. The power law decay was set to (i(v) = («;_f)~ 3 
and the functions r(v ) and f—(v) that we used are given 
by Eqs. ( 60 ) and ( 61 ) respectively, where we retained the 



first two terms in the sums. Fig. [7| displays the result of 
numerically integrating the ODE, where in the top figure 
we set the initial condition m + to 0.2 and in the bottom 
figure it was set to 0.3. The fact that these solutions 
match the DTR results is quite intriguing, and suggests 
that perhaps using Einstein's classical field equation in 
this instance is a more reliable approximation than might 
originally be expected. 
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FIG. 4: Fa plotted as a function of itia with fixed 
values r = 1, P = 0.6 and a = 0.5. The top figure gives 
a large view of the function; note that Fa diverges at 
iriA = —r/2P = —5/6. Although it isn't obvious from 
the figure, Fa asymptotes to a value of 
r 4 /(r 4 + Oq) = 0.8 when itia — > ±oo. The bottom figure 
displays the range of itia from to 3 and we include a 
horizontal line to emphasis that, given a value of Fa, 
there are two possible values of ttia ■ As Fa increases we 
see that uia will either decrease or increase depending 
on whether the initial value is taken to be on the left or 
right side of the minimum in Fa, respectively. 
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FIG. 5: Solution r(v) of an outgoing null geodesic, with 

parameters M = 1, P 2 = 0.1 and ofa = 0.15, giving 
r_ = 0.2. The solution quickly passes by the radius 0.2 
before coming to a stop (this is where it intersects the 

inner apparent horizon) and starts to increase, 
eventually settling onto the Cauchy horizon at v = oo. 
The vector field indicates the slope of r{v) at each point 
(v,r). Specifically, this is given by dr/dv = f/2. 




FIG. 6: Schematic illustrating the situation considered 
in the Ori model. We see both the continuous ingoing 
radiation as well as the outgoing shell S. This shell acts 
to split the spacetime into two regions. 
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FIG. 7: The results of numerically integrating Eq. 
with parameters M = 1, P 2 = 0.25 and a 2 , = 0.1; 
gives —r_/2P = —MP — —0.5. For the top graph the 

initial conditions were m+ = 0.2, which results in 
m+(i>), asymptotically approaching — r_/2P = —0.5. 

For the bottom graph the initial conditions were 
m + = 0.3, which results in m+(v) diverging to infinity 
at finite v. These are the two solutions obtained from 
the DTR analysis. 



